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ABSTRACT 



The J — K s color distribution (CD) with a bin size of 0.05 magnitude for the 
entire Milky Way has been carried out by using the Two Micron All Sky Survey 
Point Source Catalog (2MASS PSC). The CDs are bimodal, which has a red peak 
at 0.8 < J — K s < 0.85 and a blue peak at 0.3 < J — K s < 0.4. The colors of 
the red peak are more or less the same for the whole sky, but that of the blue 
peak depend on Galactic latitude (J — K s ~ 0.35 at low Galactic latitudes and 
0.35 < J — K s < 0.4 for other sky areas). The blue peak dominates the bimodal 
CDs at low Galactic latitudes and becomes comparable with the red peak in 
other sky regions. In order to explain the bimodal distribution and the global 
trend shown by the all sky 2MASS CDs, we assemble an empirical HR diagram, 
which is composed by observational-based near infrared HR diagrams and color 
magnitude diagrams, and incorporate a Milky Way model. In the empirical 
HR diagram, the main sequence stars turnoff of the thin disk is relatively bluer, 
(J — K s )q = 0.31, when we compare with the thick disk which is (J — K s )q = 0.39. 
The age of the thin/thick disk is roughly estimated to be around 4-5/8-9 Gyr 
according to the color-age relation of the main sequence turnoff. In general, the 
2MASS CDs can be treated as a tool to census the age of stellar population of the 
Milky Way in a statistical manner and to our knowledge this is a first attempt 
to measure the age. 

Subject headings: Galaxy: general - Galaxy: stellar content - Galaxy: structure 
- stars: luminosity function, mass function - stars: Hertzsprung-Russell and C— M 
diagrams - infrared: stars 
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Introduction 



One of the obvious properties of the Milky Way that we want to learn is its morphology. 
The fact that we reside inside the Milky Way hinders us from a global view of our own 
Galaxy which we enjoy over other galaxies. William Herschel made use of star c ount (SC ) 



to ge t a rough idea of how the Milky Way looks like in the eighteenth century (jHerschel 



17851 ). The advance in d ata collecting and 



research groups (such as 



analysis methods in late twen tieth century let 



Bahcall fc Soneira 



1980 



Gilmore &: Wyse 



19851 ) to quantify the 



major components (i.e., disks and halo) of our Galaxy by SC. They concluded that the 



pr operties of the ste 



in 



Chang et al. 



2011 



lar population in different components are distinguishable (see Table 1 
and references therein). Although many studies have been carried out, 



a common consensus on the structural parameters of the Milky Way has not been reached 



yet. One of the reasons can be at 
(i.e., limiting volumes) were used ( 


;ributed to c 


ifferent sky regions and 


imiting magnitudes 


Siegel et al. 


2002; 


Karaali et al. 


2004; 


Bilir et al. 


2006a 


b 


Karaali et al. 


2007; 


Juric et al. 


2008 


). This problem can be overcome by wide sky coverage 



SC studies which become possible in the last decade as mod ern all sky surveys b ecome 



available, such as the Two Micro n All Sky Survey (2MASS 



Skrutskie et al. 



Sloan Digital Sky Survey (SDSS; 



Yo 



Rapid Response System (Pan-Starrs 



k et al. 



Kaiser et al. 



20061 ). the 



2000 i), the Panoramic Survey Telescope & 



Explorer (WISE; 



Wright et al. 



2002j) and the Wide-field Infrared Survey 



20101 ) . However, subsequent researches did not settle the 
issue and the variation in different studies is thought to come from the degeneracy between 
the structural parameters. Moreover, the Galactic struct ural parameters have been reporte d 



to show dependence on Ga 



Cabrera-Lavers et al. 



2007 



actic longitude and latitud e ( 



Bilir et al. 



2008 



Yaz et al. 



Bi 



ir et al 



2006a 



Ak et al. 



2007 



2010). This makes the measurement 



and interpretation far more complex and challenging. 

To solve the degeneracy, information other than SC is needed, e.g., dynamics, color, etc. 
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The most handy tool is stellar color as it is available in all sky surveys. That is the reason 
why most SC studies come with color distribution (CD), which calculates the number of 
star in each color bin. In addition to luminosity function (LF) and density profile (DP) used 
in single wavelength SC study, CD requires Hertzsprung-Russell diagram (HR diagram) to 
transform each luminous star into its corresponding color. HR diagram could be empirical 
(i.e., from observation) or theoretical (i.e., from population syntheses). In the empirical 



case, HR diagram and number densiti es of each spectral 



cluster and loca 



Wainscoat et al. 



field stars (see e.g. 



Bahcall fc Soneira 



Type 



1980 



s tar are provided from ste llar 



Mamon fc Soneira 



1982 



19921 ). Difficulties include time-consuming data collection and processing, 



and possible bias of sample selection (e.g., incompleteness, improper representation). In 
general, observational HR diagram in optical wavelength is obtained first and then other 
waveband HR diagrams are constructed by color transformation (see 



Mamon &; Soneira 



1982 



Wainscoat et al 



1992 



, for examples). Thus, the uncertainty in color transformation 
needs to be considered carefully. In the theoretical case, HR diagram is deduced from 
the Hess diagram, which is generated by population synthesis based on knowledge of the 
initial mass function, ste llar evolution and star format ion history, etc. T he Besangon 



model ( Robin et al. 



20051 ) are examples. 



20031 ) and the TRILEGAL model flGirardi et al. 
Moreover, population synthesis models depend greatly on input parameters, and stellar 
evolution and atmosphere. 



In 



Chang et al 



(12010 



20111 ). we have shown that the K a all sky 2MASS SCs can be 
well described by a single power law LF and a three-component DP (namely, the thin disk, 
the thick disk and the halo). In this work, we would like to extend our model to explain 
the J — K s CDs of the entire Milky Way. We gather currently available NIR data, i.e., 
from direct NIR observations, transformations of optical HR diagram, and observational 
NIR color magnitude diagram (CMD) of star clusters, to assemble an empirical NIR HR 
diagram and incorporate a Milky Way model to explain the 2MASS data. In a way, the HR 
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diagram we used is a 2MASS optimized HR diagram. The article is organized as follows. 
The features of the all sky 2MASS CD are described in section 2. Our Milky Way model, 
including LF and DP, is provided in section 3. The procedure to create our Hess diagram 
and HR diagram is given in section 4. In section 5 we present our results on CDs of the 
Milky Way, and section 6 is a summary and concluding remarks. 



2. The 2MASS Data and Its Color Distributions 



We use the 2MASS Point Source Catalog (2MASS PSC, ICutri et all 120031 ) to carry 
out J — K s CDs with a bin size of 0.05 magnitude for the entire Milky Way. We select 
objects with the following criteria: (1) signal-to-noise ratio > 5, (2) detection in all J, 
H, K s bands and (3) K s magnitude between 5 and 14 mag. The last criterion ensures 
99% completeness rate before the K s limiting magnitude (i.e., 14.3 mag) and avoids the 
relatively large photometric error for K s < 5 mag objects. The whol e sky is divided into 



8192 nodes according to level 5 Hierarchical Triangular Mesh (HTM, iKunszt et al.l 120011 ). 
which samples the whole sky roughly evenly and has a 2 degrees angular separation on 
average. The radius of each node is 1 degree (i.e., each node covers ir square degree). 
Because of the shallower limiting magnitudes and the complex extinctions in the Galactic 
center area, this work does not attempt to explain the 2MASS CDs in this region. 

Several common features are identified in the whole sky 2MASS CDs, which are listed 
below and demonstrated in Fig. [TJ 



1. Most of the 2MASS CDs are bimodal, which has a blue peak at 0.3 < J — K s < 0.4 
and a red peak at 0.8 < J — K s < 0.85. The 2MASS CDs at low Galactic latitudes 
have one more peak at 0.55 < J — K s < 0.6. Fig. [TK shows three examples of typical 
2MASS CD at different Galactic latitudes. 
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2. The blue peak dominates the 2MASS CD at low Galactic latitudes and becomes 
comparable to the red peak at high Galactic latitudes. The middle peak only shows 
up at low Galactic latitudes and does not exist in other sky areas any more. Fig. [Tb 
shows the averaged CDs over 225° < I < 255° at different Galactic latitudes. 

3. The J — K s colors of the all sky red peaks and the middle peaks at low Galactic 
latitudes are almost fixed. However, the colors of the blue peak depend slightly 
on Galactic latitudes, which is J — K s ~ 0.35 at low Galactic latitudes and 
0.35 < J — K s < 0.4 at medium and high Galactic latitudes. The situation is 
demonstrated in Fig. [Tb, which is the same as Fig. [lb except that the total count is 
normalized to b = 10°. In contrast, the color of the blue peak does not change along 
Galactic longitude, see Fig. [TH & [Tb. 

4. We call the tapering off distribution on the blue side of the blue peak the blue wing. 
The shapes of the blue wings at different Galactic latitudes are very similar. Due to 
their bluer blue peaks, the blue wings of the CD at low Galactic latitudes can extend 
to J — K s = 0. However, the blue wings at medium and high Galactic latitudes 
seldom extend to less than J — K s = 0.2. The normalized 2MASS CDs in Fig. [Tb 
clearly shows this property. 

5. Occasionally, some 2MASS CDs at medium or high Galactic latitudes have a middle 
peak, which is different from the middle peak at low Galactic latitudes. We call it 
'the extra peak' and will explaine more in section 14.11 and 14.21 The solid lines in the 
first column of Fig. [2] show three examples of the extra peak. 

6. The shape of the 2MASS CD can be altered significantly by severe extinction, which 
is the usual case at low Galactic latitudes. Fig. [3] shows the 2MASS CDs at |6| ~ 16°. 
The figures are arranged according to their E(J — K s ) extinction values. As the 
E(J — K s ) values increase, the sharpness of the 2MASS CDs decrease gradually, and 



-7- 



the J — K s colors of the blue and middle peaks become redder and redder. Moreover, 
the tail in the red end is elongated. 



3. The Milky Way Model 



Our Milky Way mo del has a three-component DP and a single power law LF 



(jChang et al. 



2010 



201ll ). The three-component DP n(R, Z) includes a thin disk D\, a 



thick disk D 2 and an oblate halo S, 



n{R, Z) = n [D^R, Z) + D 2 {R, Z) + S(R, Z)\ , ( 

where R is the galactocentric distance on the Galactic plane, Z is the distance from 
the Galactic mid-plane and no is the local stellar density of the thin disk at the solar 
neighborhood. 

The disks are in a double exponential decay form in which stellar density decreases 
exponentially along R and Z, 

(R — Rq) (\Z\-\Z e \) 



Di(R, Z) = /-exp 



(2) 



Hri Hzi 

where (Rq,Z q ) is the location of the Sun, H ri is the scale-length, H zi is the scale-height, 
and fi is the density ratio to the thin disk at the solar neighborhood. The subscript % — 1 
stands for the thin disk and % = 2 stands for the thick disk. 



The halo is a power law decay oblate spheroid flattening in the Z direction, 

W) = /J* 2 + (Z/ * )21 " 2 



(3) 



R% + {Z q /kY 

where k is the axis ratio, p is the power index and f\ is the local halo-to-thin disk density 
ratio. 
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The Mk s LF is a single power law ( jChang et al.ll2010l ). 

2icg„io (7 -i) 



10 2( 7 -l)A//5 



5 [lO^f-^Mf/ 5 - 102(7-l)M b /5j 

where 7 is the power law index, M& and Mf are the bright and faint cutoffs, respectively. 



(4) 



Note that ip(M) includes all luminosity classes, 
of the bright end and f aint end distributio ns of 



'or convenience, we adopt the peak values 



Mf = 6.5 (see Fig. 6 in IChang et al. 
note that the averages in 



Chang et al 



Chang et al. 



(I2010h . i.e., M b = -8 and 



20101). as the bright and faint ends in this study. We 
J2OIO) are M b = -7.86 ± 0.60 and M f = 6.88 ± 0.66. 



Nevertheless, the model prediction shows no significant difference between using the peak 
values and the average ones. 



The extinction model is adopted from the new COBE/IRAS result 



and the color excess is E(J — K S )/E(B 



parameters taken from 



Chang et al. 



told) and 



V) = 0.53 flSchlegel et a 



Chang et al. 



Chen et al 



1999|) 



1998|). Tabled] lists the 



(1201 lh (hereafter MW-I). 



4. The Empirical NIR HR Diagram Optimized for the 2MASS CD 

We describe the search for an empirical NIR HR diagram (tailored for 2MASS data) in 
three subsections for a simple and clear exposition. 



4.1. Basic NIR HR Diagram 

We take observational-based NIR HR diagrams in the literature as our first construction 
of the loci of ma i n sequ ence (MS) and giant branch (GB ) in M k-{J — K s )o HR diagram. 



Wainscoat et al. 



( 1992 



hereafter W92) and 



Covey et al. 



(120071 . hereafter C07) have the 



information of NIR HR diagram (i.e., the NIR color, absolute magnitude and spectral 
type). W92 also gives the number density in the solar neighborhood and the normalization 
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in different components of Milky Way. Fig. H] shows the data points of MS, GB and 
super giants taken from W92 and C07. Both data sets agree to each other very well and 
shows several common features: (1) a sharp downturn around the end of MS (hereafter, 
MS-downturn), (2) a slight upturn around late A type star (i.e., Mk s ~ 3 mag and hereafter, 
MS-turnoff), and (3) a sharp curve on GB around Mk s ~ — 1 mag (hereafter, GB-curve). 
We fit W92 data to find the loci and integrate it with MW-I to obtain a synthetic CD. 
When the synthetic CD is compared with 2MASS CD, some significant discrepancy is 
observed: (1) over predictions for the middle peak and for the blue wing (i.e., J — K s < 0.3 
mag), (2) under prediction on the blue peak, and (3) the locations of the blue and middle 
peaks are not correct. The first column in Fig. [5] demonstrates several results for different 
Galactic latitudes. However, we obtain a very useful information that the J — K s colors 
of each peak correspond to that of each sharp turn (or 'vertical' sections) on the loci. 
Apparently, the blue, red and middle peaks are the results of the MS-turnoff around late A 
type stars, the MS-downturn and the GB-curve, respectively. 



4.2. The Optimization for the Empirical NIR HR Diagram 



To study the sharp tur ns in more deta il , we u t ilize several NIR CMD s of star cluste rs 



from the literature, such as 



Beletsky et al. 



fl2009[ ): 



Saraiedini et a 


. ( 


200< 


)); 


Troisi et al. 


Beletskv et al. 


2009 




Saraiedini et al. 



1120101). Some of them (e.g., M67 and NGC6791 in 
20091 ). which have information down to the faint end of MS, are used to conclude the 
important features for our analysis: (1) the MS-downturn (J — K s )q ~ 0.83 is a common 
feature at the faint end of MS; (2) the MS-turnoff (i.e., the slight upturn around late A 
star) reflects stars departing from the MS at the end of their MS stage; (3) the MS-turnoff 
always accompanies a AM^ s ~ 2.5 K s mag vertical height extension and its (J — K s ) 
color depends on the age of the star cluster; and (4) the main sequence belt between the 
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MS-turnoff and the MS-downturn is approximately linear. Since only few GB stars are 
shown on these NIR CMDs, we are not able to get more details of the GB-curve. However, 
we believe that the location of the GB-curve should be similar to that of W92 and C07. 

Accordingly, a working HR diagram includes the following features. The MS is divided 
into three parts: (1) early MS (i.e., O-A type stars and Mk„ < 3 mag); (2) nearly linear 
late MS (i.e., F-K type stars and 3 < Mk s < 5 mag); and (3) sharp and almost vertical 
MS-downturn at (J — K s ) ~ 0.83 (i.e., M type and later stars and Mx a > 5 mag). The 
GB is separated into up and low GBs. Both are linear and connected by the GB-curve 
at (J — K s ) ~ 0.6 and M^ s ~ — 1. Consequently, we create an empirical M^ a against 
(J — K s )o locus including a vertical MS-downturn, a linear late MS, a MS-turnoff with 
a AMx a = 2.5 K s mag vertical extension, a linear low GB, a linear up GB and a flat 
transition between the tip of the MS-turnoff and the low GB. The slopes of each linear 
part (i.e., the up GB, the low GB and the late MS) and the connecting joints of different 
parts are taken from the CMDs mentioned above and the data in W92 and C07. The 
early MS is not included because the 2MASS CDs are generally deficient in number of 
stars in the region (J — K s )o < 0.2. Since the J — K s colors of the red peak of 2MASS 
CD is between 0.8 to 0.85 and that of the MS-downturns of M67 and NGC6791 are at 



0.83 f lBeletskv et al. 



2009 



Sarajedini et al. 



20091 ). we set (J - K 8 ) = 0.83 for the 



vertical MS-downturn. Moreover, the middle peak is around J — K s = 0.55 to 0.6 and 
the GB-curves in W92 and C07 also have similar color, so we set (J — K s )q = 0.6 for the 
GB-curve and let the low GB between (J — K s )o = 0.5 and 0.6. To incorporate the idea 
that different (J — K. q ) n colors of the MS-turnoff for different evolutionary stage populations 



( jSarajedini et al. 



20091 ). we allow the MS-turnoff to move along the main sequence belt. 



The left panel of Fig. H] shows the empiric al NIR locus on top of the data points taken from 



W92, C07 and M67 f lBeletskv et al 



20091 ). To mimic observational color scattering, we put 



a 0.06 color dispersion to the empirical locus (see the right panel of Fig. H]). 
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Incorporating this empirical NIR HR diagram with MW-I, we found several features 
in the model prediction (see Fig. E]): (1) the CD generated by the thin disk is a bimodal 
distribution; (2) the CD generated by the thick disk only shows the blue peak; (3) both CDs 
of the thin and thick disk can have a middle peak at low Galactic latitudes; (4) the overall 
model prediction is mainly dominated by the CD of the thin disk; (5) the blue peak of the 
thick disk has an appreciable contribution to that of the overall model prediction only at 
medium and high Galactic latitudes; and (6) the red peak of the overall model prediction is 
determined mostly by the thin disk. 

At first, we focus on medium and high Galactic latitudes to avoid severe extinction. 
The (J — K s )o colors of the MS-turnoff for the thin and thick disk are set between 0.35 
and 0.4 to comply with the blue peak of 2MASS CD at medium and high latitudes. We 
integrate the empirical locus (Fig. H]) with our Milky Way model MW-I to generate the 
synthetic CD. For reasonable good fit, the (J — K s ) color of the thin disk and the thick 
disk are found to be 0.34-0.35 and 0.38-0.39, respectively. The second column in Fig. [5] 
shows some results at different Galactic latitudes. 

The result at medium and high Galactic latitudes agrees well with the 2MASS data. 
Then, we inspect low Galactic latitudes which have relatively small extinction. The result 
shows significant difference in the blue wing of the blue peak (i.e., J — K s < 0.3), but the 
rest of the CD is still in a good agreement (see the first row of the second column in Fig. [5]). 
If a bluer turnoff is assigned to the thin disk, the discrepancy in the blue wing can be 
compensated somewhat, but it will cause an over prediction in the blue wing and a shift in 
J — K s color of the blue peak at high Galactic latitudes (see the third column in Fig. [5]). 
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4.3. The Final Adjustment 



When experimenting parameters of the model in previous subsections, we notice a 
small scale height thin disk would make its blue peak change more dramatically along 
Galactic latitude. This makes the overall blue peak (i.e., the thin and thick disks ) at low 
Galactic latitudes be dominated by the thin disk, while the blue peak at medium and high 
Galactic latitudes be dominated by the thick disk. Therefore, a small scale height thin 
disk with a bluer MS-turnoff not only can eliminate the difference in the blue wing at low 
Galactic latitudes, but also can reduce its contribution to the blue peak at medium and 
high Galactic latitudes. In this way, the J — K s color of the blue peak could be relatively 
bluer at low Galactic latitudes and relatively redder at the other sky region. With this 
strategy, we select an acceptable Galactic structure co nfiguration with a r elatively smaller 



scale height thin disk from our previous SC study (see 



Chang et al. 



2011 



for details) and 

come up with the best fit model as follows: (1) the scale height of the thin disk is H z i = 260 
pc and other structural parameters are listed in Tabled] (hereafter, MW-II); and (2) the 
MS-turnoff of the thin disk is (J - K s ) = 0.31 mag. With this model, the 2 MASS CDs at 
low Galactic latitudes can be fitted much better than that in subsection 14.21 (see the last 
column in Fig. |5]). 

Fig. |4] shows the HR diagram obtained fro m subsections 14.21 and 14.31 Incorporating the 



luminosity function derived from 2MASS SC ( jChang et al. 



2010|), we construct the Hess 



diagram of MW-II in the right panel of Fig. |H The joints on the HR diagram of MW-II are 
listed in Table [2j 
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5. Results and Discussion 

In Figs. [6]&[7]we compare our model prediction (i.e., MW-II in section H~3j) and the 
2MASS CDs at different Galactic latitudes. In order to minimize the uncertainty caused by 
extinction correction at low Galactic latitudes, we compare sky areas with relatively small 
extinction. The model prediction and the 2MASS CD agree very well in the overall shape 
and the J — K s color of each peak. The number variations along Galactic latitudes of each 
peak are well reproduced by the model. 

In order to present the global trend of the 2MASS CDs, we show the number ratio of 
the blue part (i.e., J — K s < 0.6) to the red part (i.e., J — K s > 0.6) in Fig. [8] along with 
the prediction of MW-II. Both figures have very similar trend which is relatively blue at 
low Galactic latitudes and become comparable in other sky regions. The very red Galactic 
disk is the result of severe extinction. 

When we inspect the Galactic disk region (i.e., |6| < 10°), the model prediction has 
obvious discrepancy with the 2MASS CD. The discrepancy can mostly be attributed to 
improper extinction correction and number inconsistency between the model prediction 



and t 



re 2MASS data. The number inconsistency, which has been reported in 



Chang et al. 



( 120 1 lh . results in gaps on each peak between the model prediction an d the 2MASS CD. 



Chang et al. 



(1201 ih . 



Fig. |9] demonstrates some cases together with the SC results taken from 
In these cases, the sky areas have relatively small extinction. This helps us to delineate the 
discrepancy mainly comes from the number inconsistency (see the difference in SC result, 
lower panels in Fig. [9]). Several examples of improper extinction correction are illustrated 
in the first column of Fig. [TUJ The improper extinction can make a significant offset on the 
J — K s colors of the blue and middle peaks and distort the shape of CD. Moreover, we also 
find a significant blue population (i.e., J — K s < 0.2) in the 2MASS CD on the Galactic 
disk < 10°) around the Galactic longitudes of 180° < / < 240° (see the second column 
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of Fig. [TO]) . This blue population is neither predicted by the model nor the blue wing of the 
blue peak. If we omitted the extinction correction, the model still cannot account for such 
a blue population. The origin of this blue population warrants further studies. 

Several 2MASS CDs at medium and high Galactic latitudes have an extra peak at 
J — K s ~ 0.55. Their occurrence and height do not correlate with Galactic coordinates. 
Most of them are insignificant, but we show three good examples in the first column of 
Fig. |2] (solid line). We note that each example has a globular cluster within the field of 
view. In order to nail down the origin of the extra peak, we generate globular-cluster-only 
CD (dashed line). Each globular-cluster-only CD has a spike, which has J — K s color very 
similar to that of the extra peak in the whole field CD (i.e., the CD contains stars both 
from the globular cluster and the field; solid line). We also obtain the CD of the ambient 
field (dot-dashed line) by subtracting the globular-cluster-only CD from the whole field 
CD. The extra peak does not show up anymore in the ambient field CD. Besides, each 
globular cluster CMD has an obvious GB with an obscure GB-curve at J — K s ~ 0.55 
similar to that on the empirical HR diagram (see Fig. H]). Therefore, we conclude that the 
extra peak is due to the low GB of the globular cluster and its significance depends on 
the relative amount of star in the low GB of the globular cluster to that in the ambient 
field. Besides, if the limiting magnitude could be reached down to the MS-turnoff of the 
globular clusters, the observed blue peak would have more stars than the model prediction. 
Similar number enhancement could happen in the 2MASS CDs if there is an overdensity 
within the field of view, which have a dominant stellar population over the smooth stellar 
distribution of the Milky Way (e.g., star cluster, stellar stream, Galactic bulge, Galactic 
arm, etc.). For instance, the extremely dense Galactic bulge, which is not included in our 
model, makes a prominent middle peak on the 2MASS CDs, which obviously outnumbers 
the model prediction (see the third column of Fig. [TOl) . However, we do not observe any 
obvious number enhancement in other 2MASS CDs which have known stellar streams in 
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their field of view, e.g., 



( jPenarrubia et al 



2005 



Sagittarius dwarf galaxy ( IMajewski et al.l 120031 ) . Monoceros stream 



Table 1), etc. It is possible that within the 2MASS detection 



limit the populations of these overdensities are too small. 

At low Galactic la titudes there is another possible way to generate the middle peak. It 



is the red clump stars ( lAlves 



2000 



Cabrera-Lavers et al. 



20071 ) . which occupies the space 



close to the low GB on our HR diagram and has a higher number density relative to the 
other stars on GB. Such relatively higher number density would produces more stars in 
the corresponding color bins. In our model we do not include the red clump stars in our 
LF and HR diagram, but our model still have good agreement with the 2MASS CD. The 
reason might be the net contribution of the red clump stars in generating the middle peak 
is similar to that of the low GB. 



Our best model that can explain the all sky 2MASS CDs (see section I4.3P comprises 
a thin disk, which has a relative short scale height (260 pc) and a relative blue MS-turnoff 
(i.e., younger), and a thick disk (1040 pc), which has a relative red MS-turnoff (i.e., 
older). Our best model, MW-II, that can explain the all sky 2MASS CDs (see section FOj) 
comprises a thin disk (260 pc) with a relatively blue MS-turnoff (0.31 mag, i.e., younger), 
and a thick disk (1040 pc) with a relatively red MS-turnoff (0.38 mag, i.e., older). The 
two populations can be interpreted in the context of the formation of the Milky Way. The 
thick disk was formed in an epoch earlier than the thin disk. Moreover, both disks have 
the same single power law LF which implies that the initial mass function does not change 
along the Milky Way's evolution (unless the mass-luminosity relation was time dependent 
in the past). In this sense, the all sky 2MASS CDs can be used as a tool to census the age 
of stellar popula t ion in the Milky Way. If we apply the color-age relation of MS-turnoff of 



Sarajedini et al. 



( 120091 ) to our empirical HR diagrams, we can roughly estimate the ages 



of the thin and thick disk to be 4-5 Gyr and 8-9 Gyr, respectively. The age estimation 
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should be viewed as the lower limit. The blue boundary of the 2MASS CDs do not allow 
bluer MS-turnoffs of the disks. Moreover, the 0.05 mag bin size of our 2MASS CD can only 
provide a low color resolution of the blue peak, hence the color of the MS-turnoff and the 
derived age cannot be estimated more precisely. Nevertheless, to the best of our knowledge 
our result could be the first attempt to measure the age of global stellar components in a 
general statistical sense. We conclude that the whole population of the thick disk is about 
3-4 Gyr older than that of the thin disk. 

In this work, we seldom mention the halo component because it has very 
limited contribution to the total CDs in magnitude range we are interested (i.e., 
4 mag < K s < 14 mag). Thus, we are not able to extract solid and useful information of 
the halo from 2MASS CD. 

6. Summary 

We use 2MASS PSC to carry out the whole sky J — K s CDs, which show a bimodal 
distribution with a blue peak and a red peak. The blue peak is the dominant feature at 
low Galactic latitudes and becomes comparable to the red peak at high Galactic latitudes. 
The all sky J — K s colors of the red peak are almost the same (between 0.8 and 0.85 mag). 
However, the J — K s colors of the blue peak at low Galactic latitudes are a bit bluer (0.35 
mag) than that at medium and high Galactic latitudes (between 0.35 mag and 0.4 mag). 
However, the J — K s colors of the blue peak at low Galactic latitudes are a bit bluer (which 
is J — K s ~ 0.35) than that at medium and high Galactic latitudes (which is between 0.35 
and 0.4). Besides, a middle peak shows up at J — K s ~ 0.55 at low Galactic latitudes and 
does not exist in other sky areas. Several medium and high Galactic latitudes have an extra 
peak at J — K s ~ 0.55, which is not the same as the middle peak at low Galactic latitudes. 
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In order to explain the whole sky 2MASS CDs, we create an empirical Mk s against 



J-K s HRc 
and CMDs ( 


iagram optimized for t 


he 2 MASS data by 


Wainscoat et al. 


1992: 


Covev et al. 


2007: 


2009; 


Troisi et al. 


2010 


) and incorporate a Milky Way 



gathering: avai 


able 


sTIR HR diagrams 


Beletskv et al. 


2009 


Sara 


ledini et al. 


model f 


Chang et al. 


2010. 


2011 


). The 



HR diagram of the thin disk has a bluer MS-turnoff (i.e., younger), and that of the thick 
disk has a redder MS-turnoff (i.e., older). For a better explanation of the 2MASS CDs, the 
result is in favor of a relatively small scale height thin disk (260 pc, MW-II, compare with 
360 pc in MW-I; other parameters of MW-I and MW-II are given in Tabled]). We note 
that the configuration of MW-II is within the 'acceptable' solutions due to the degeneracy 



between the structural parameters, w : 
to study the Galactic structure (e.g. 



rich has been poin ted out as a problem for using SC 



Chang et al. 



201ll ). Now information from CD helps 



to rule out the degeneracy, and we deem that MW-II should be the preferred model. 

We find that the blue peak, the red peak and the middle peak are due to the MS-turnoff, 
the MS-downturn and the GB-curve, respectively (see Fig. H]). Moreover, the extra peak at 
some medium or high Galactic latitudes is due to globular clusters. Our model cannot trace 
the CDs which suffer severe extinction and need a more elaborate extinction correction. 



When we apply the color- age relation of the MS-turnoff (jSarajedini et al.ll2009l ) to our 
empirical HR diagram, the thin and thick disks are estimated to be around 4-5 Gyr and 
8-9 Gyr, respectively. The idea is consistent with the Milky Way's formation theory. We 
realize that the whole sky 2MASS CDs can be used as a tool to census the age of stellar 
populations of the Milky Way in a statistical manner. To the best of our knowledge, this 
study is the first attempt to measure the ages of the thin and thick disks from the whole 
sky stellar population. 

It is no doubt that the 2MASS PSC provides an unique tool to study the global 
property of the Milky Way, not only for its tremendous sky coverage (which avoids the 
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selection effect of limited sky coverage), but also for the benefits of using NIR wavelength 
(which keeps fine angular resolution and is less affected by interstellar extinction). 

We acknowledge the use of the Two Micron All Sky Survey Point Source Catalog 
(2MASS PSC). This work is supported in part by the National Science Council of Taiwan 
under the grants NSC-98-2923-M-008-001-MY3 and NSC-99-2112-M-008-015-MY3. 
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Table 1. The Milky Way model. 



MW-I 



MW-II 



Density Profile 
Thin Disk 

H r l 

H zl 

no 

Z @ 

H: 

Thick Disk 
H r2 

h 

Spheroid 
K 
P 
fh 

Mk s Luminosity Function 

7 

M b 
Mf 

Extinction Correction 

new COBE/IRAS result (Chen et al. 1999) 
E(J — K 3 )/E(B — V) 



3.7 kpc 
360 pc 

0.030 stars/pc 3 
25 pc 

8 kpc (Reid 1993) 

5.0 kpc 
1020 pc 
7% 

0.55 
2.6 

0.20% 

1.85 

-8 

6.5 



0.53 (Schlegel et al. 1998) 



3.7 kpc 
260 pc 

0.039 stars/pc 3 



5.0 kpc 
1040 pc 
10 % 
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Table 2. The Empirical HR Diagram Optimized for 2MASS. 



Part (J - K 3 ) M Kb 



Common Joints 

Faint End of MS 0.83 6.5 

MS-Downturn 0.83 5.0 

GB-Curve 0.6 -1.0 

Bright End of Up GB 1.4 -8.0 

Thin Disk 

MS-Turnoff 0.31 3.34 

Bright End of MS-Turnoff 0.31 0.84 

Faint End of Low GB 0.54 0.34 

Thick Disk 

MS-Turnoff 0.38 3.56 

Bright End of MS-Turnoff 0.38 1.06 

Faint End of Low GB 0.53 0.56 
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Fig. 1.— Salient features on the 2MASS CD. (a) Typical examples of 2MASS CD, [I, b] = 
[235.6°, 8.1°] (solid line), [248.5°, 20.1°] (dashed line) and [243.0°, 57.3°] (dotted-dashed line), 
(b) The 2MASS CDs averaged over 225° < I < 255° at different Galactic latitudes. The 
Galactic latitudes are indicated by different colors, (c) Same as (b), except that the total 
count is normalized to that of b = 10°. (d) The 2MASS CDs averaged over 15° < b < 20° 
at different Galactic longitudes. The Galactic latitudes are indicated by different colors, (e) 
Same as (d), except that the total count is normalized to that of / = 60°. 
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Fig. 2. — The first column is the 2MASS CDs having the extra peak at medium and high 
Galactic latitudes. The whole field, the globular-cluster-only and the ambient field CDs are 
represented by the solid, dashed and dotted-dashed lines, respectively. The identification 
of the star cluster and its center coordinate are indicated on the upper-right corner of each 
figure. The second, third and last columns show the CMDs of the whole field, the globular- 
cluster-only and the ambient field, respectively. The radii of the field of view of the whole 
field and the ambient field CMDs are 1 degree and that of globular-cluster-only CMD is 
given on the upper-right corner of each figure in the third column. 
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Fig. 3. — The 2MASS CDs around |6| ~ 16° with different E(J — K s ) extinction values. The 
legend gives the Galactic coordinates and the E(J — K s ) value. 
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Fig. 4. — The HR diagram and the Hess diagram. The left panel is the NIR HR diagram. 
The data points are taken from W92 (color symbols), C07 (black symbols) and M67 (small 



black open circle; 



Beletsky et al 



2009I ). The upper- right legend indicates luminosity classes 
and the number density is coded by the color bar. C07 does not have number density 
information, so it is plotted as black. Each small black open circle represents one member 
star of M67. The lines represent different HR diagrams used in subsections I4.2[ I4.3I The 
upper-left legend gives the corresponding models and components. The right panel is the 



Hess diagram of the thin disk used in MW-II whose number density is color coded. 




Fig. 5. — The synthetic CDs of different HR diagrams (see section H] for details). Black line 
is the 2MASS data, and pink line is the model prediction. The blue, green and orange dotted 
lines represent the thin disk, the thick disk and the halo, respectively. 
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Fig. 6. — Comparison between 2MASS CDs (black) and model predictions (pink). The 
coordinates and the extinction E(J — K s ) are given at the upper-right corner of each figure. 
The blue, green and orange dotted lines represent the thin disk, the thick disk and the halo, 
respectively. 
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Fig. 7. — Continuation of Fig. [6j 




Fig. 8. — The number ratio of the blue part (J — K s < 0.6) to the red part (J — K s > 0.6). 
Top: 2MASS data. Bottom: model prediction. Note that the very red Galactic disk is the 
result of severe extinction and it does not reflect the intrinsic number ratio in these regions. 
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Fig. 9. — The number inconsistency between the model prediction and the 2MASS CD. The 
upper panel is the comparison of CDs. The upper-right corner shows the coordinates and 
the extinction E(J — K s ) value. The solid-black, solid-pink, dotted-blue, dotted-green and 
dotted-orange represent the 2MASS CD, the model prediction, the thin disk, the thick disk 
and the halo, respectively. The lower panel is the comparison of star counts between 2MASS 
data (plus) and the model prediction (line). 
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Fig. 10. — The first column shows the improper extinction correction. The second column 
is the examples of blue population. The third column is the examples in Galactic center re- 
gion. The solid-black, solid-pink, dotted-pink, dotted-blue, dotted-green and dotted-orange 
represent the 2MASS CD, the model prediction, the model prediction without extinction 
correction, the thin disk prediction, the thick disk prediction and the halo prediction, re- 
spectively. 



